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ABSTRACT 

We show that the observed relationship between the fraction of low-mass X-ray binaries (LMXBs) 
found in globular clusters (GCs) and the GC-specific frequency for early-type galaxies is consistent with 
an LMXB formation model in which the field population of LMXBs is formed in situ via primordial 
binary formation. The suggestion that a significant fraction of the field LMXB population in early- 
type galaxies was formed in GCs is not required by the data. Finally, we discuss observational studies 
that will test this model more thoroughly. 

Subject headings: binaries: close — X-rays: binaries — X-rays: galaxies — galaxies: elliptical and 
lenticular — globular clusters: general 



1. INTRODUCTION 

Low angular resolution studies of early-type galaxies 
found that the X-ray emission from these systems typi- 
cally consists of two components, a soft, thermal com- 
ponent from hot gas, and a harder, power-law com- 
ponent attributed to unresolved X-ray binaries (e.g., 
Matsumoto ct al. 1997). The superb angular resolution 
of the Chandra X-ray Observatory has allowed for a sig- 
nificant fraction of the power-law component to be re- 
solved into point sources, c onfirming the binary nature 
of the hard component (e.g.. ISarazin et al]l2000li . Early- 
type galaxies are generally dominated by stellar popu- 
lations older than 1 Gyr. Thus, the majority of X-ray 
binaries in these galaxies are low-mass X-ray binaries 
(LMXBs), rather than the short-lived (~ 10 7 yr) high- 
mass variety. LMXBs consist of a neutron star (NS) or 
black hole (BH) accreting material from a < 1 M Q donor. 

A number of early-type galaxies have been studied with 
Chandra in an effort to understand the LMXB popula- 
tion of these galaxies. Combined X-ray and optical stud- 
ies have found that 20%-70% of the LMXBs in early- 
type galaxies are associated with globular clusters (GCs 
Angelini et al.l 120011 iBl anton et all l2001t iSarazin et al 
20011 120031 iKundu et all 120021 120031 IMaccarone et all 



20031 ISivakoff et a I 120031 f IHurnnhrey fc Bnotd 12004 
Jordan et all 12004 iRandall et all 12004ft . This is signif- 



icantly larger than the roug hly 10% of LMXB s found 
in GCs in the Milky Way. IWhite et all <|2002f) . using 
ASCA data, found that the X-ray to optical luminosity 
ratio, Lx/L opt , was proportional to the GC-specific fre- 
quency, S n (the number of GCs per unit optical luminos- 
ity), of the galaxies. They found no correlati on between 
the st ellar age of the galaxies and Lx/L opt . IWhite et al.l 
( 2002) concluded that most or all of the LMXBs in early- 
type galaxies were formed in GCs and further suggested 
that the LMXBs not associated with GCs either had 
been ejected from GCs or were formed in GCs that had 
since disrupted. The similarities in the properties of the 
GC LMXBs and those found in the fields of early-type 
galaxies has been used as support for this scenario (e.g., 
IMaccarone et al.ll2003|) . 

However, it has been noted that the fraction of 
LMXBs found in GCs shows some dependence on Sn 



ijMaccarone et alJl2003i ISarazin et aDl2003|) . In this Let- 
ter, we show that the relationship between the fraction 
of LMXBs found in GCs and Sn is consistent with the 
field population of LMXBs being formed via primordial 
binary formation in the fields of the galaxies. 

2. MODEL OF LMXB FORMATION 

We assume that LMXB formation in early-type galax- 
ies is the same as commonly accepted for the Milky Way. 
In this scenario, the field LMXBs are formed in situ via 
primordial binary formation and the LMXBs associated 
with GCs are formed through dynamical interactions in 
the clusters. For any galaxy, the fraction of LMXBs 
found in GCs, denoted by /lmxb.gc, is equal to 
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Almxb.gc + -Ylmxb.f ' 



(1) 



where A^lmxb.gc an d -ZVlmxb,f are the numbers of 
LMXBs found in GCs and the field, respectively. 

It has long been recognized that GCs produce ~ 100 
times more LMXBs per unit mass than the field of the 
Galaxy. This enhanced formation efficiency is attributed 
to dynamical interactions between NSs and st ars in the 
dense environments of GCs (e.g.. lCla"rklll975j) . Obser- 
vational evidence suggests that « 4% of GCs have a 
bright {L x > 5 x 10 37 erg s" 1 ) LMXB regardless of 
the h ost galaxy morphological type (e.g., IKundu et alJ 
2002). Given the simila rities in the average pr operties of 
GCs in all galaxies (e.g.. lAshman Sz Zepll998|) . this con- 
stant LMXB formation efficiency is not surprising and is 
a measure of the efficiency of dynamical formation mech- 
anisms. We then expect that -/Vlmxb,gc is proportional 
to the number of GCs in the galaxy. The number of GCs 
is equal to Sn times the galaxy luminosity in the appro- 
priate units. For our purposes, it is more interesting to 
give -ZVlmxb.gc m terms of the galaxy mass, so we as- 
sume a fixed mass-to-light ratio for our sample galaxies. 
Therefore, we find that 



A^lmxb.gc oc SnM+ , 



(2) 



where M* is the stellar mass of the galaxy. Here, and 
throughout, we neglect differences in the stellar mass- 
to-light ratio of our sample galaxies as we are primarily 
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concerned with early-type systems. Were we to compare 
early- and late- type galaxies, it would be more appro- 
priat e to use the mass norm alized GC-specific frequency 
Ce.g.. lAshman fc Zepflll998|) . 

LMXBs in the field of the Galaxy are the descendants 
of primordial binaries containing a massive star. LMXBs 
have a long (~ 1 Gyr) lag tim e between star forma- 
tion and the X-ray active phase ijWhite fc Ghoshlfl998f) 
and therefore track the star formation history over long 
timescales. For a single starburst model for star forma- 
tion in early- type galaxies, the number of LMXBs formed 
from primordial binaries will be a function of the stellar 
mass involved in the star formation event and the time 
since star formation occurred. Initially, the number of 
active LMXBs will increas e until ~ 1 Gyr and th en be- 
gin to decrease with time ijWhite fc Ghoshl fl998'). The 
exact time dependence of the LMXB population is un- 
known. In addition, complex star-formation histories, 
such as small starbursts associated with mergers, will 
complicate the time dependence for any galaxy. For 
our model, we assume that the time dependent varia- 
tion of the field LMXB population is small when compar- 
ing different early-type galaxies. Under this assumption, 
ATlmxb,f is simply proportional to M+. 

Plugging the relationships for 7V LM xb,gc and Almxb,f 
into Equation ^ and rearranging, we find 

/lmxb,gc = f 1 + ' (3) 

where C is a constant that accounts for the formation ef- 
ficiencies in GCs and the field, and the mass-to- light ratio 
of the galaxies. To test this model, we have accumulated 
data from the literature for nine early-type galaxies with 
both Chandra observations and GC identifications. Ta- 
ble Ogives the values of Sn and the calculated fraction 
of LMXBs found in GCs. We then fit these data with a 
function of the form 

/lmxb,gc = (a + , (4) 

where A and B were both allowed to vary. We find 
that A = 1.19 ± 0.19 and B = 3.9 ± 1.0, with a re- 
duced \ 2 for the fit of 0.74. The data and best-fit func- 
tion are shown in Figure ^ For reference, we have in- 
cluded a data point for the Milky Way with S/y = 0.5 
and /lmxb.gc ~ 0.1 l|Harrislll99lULiu et alJl200ll) . The 
Milky Way data point includes both the disk and bulge 
contributions. As we see, the data follow the expected 
trend from our model. Therefore, we conclude that the 
field population of LMXBs in early-type galaxies was 
formed in the field from primordial binaries. While it 
is intriguing that the best-fit model also fits the Milky 
Way data point, it should be noted that for a luminosity 
limit similar to those found in studies of early-type galax- 
ies, the fraction of LMXBs in GCs in the Milky Way is 
highly uncertain because of low number statistics, diffi- 
culties in determining the distance to Galactic LMXBs, 
and the issue of including transient systems in the statis- 
tics. As noted earlier, mass-to-light ratio differences may 
also be important when comparing the GC-specific fre- 
quency across morphological types. 

We can also compare the prediction of this model to 
other observed properties of the LMXBs in early-type 
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Fig. 1. — Fraction of the LMXB population found in GCs (given 
in percent) plotted versus the GC-specific frequency, Sjv for 9 early- 
type galaxies. The best-fit model is overplotted. Included is a data 
point for the Milky Way LMXB population which was not used in 
the fit. 



galaxies. W hite et al.1 l|2002|) found that Lx/L opt for 
the hard spectral component in early-type galaxies was 
proportional to Sjf A . Assuming that the total X- 
ray luminosity of the LMXBs can be approximated by 
ix.total = (^x)A r LMXB, where (L x ) is an average X- 
ray luminosity per binary, assuming the same average 
lumin osity for GC and field systems fsee iSarazin et alJ 
2003), and A^lmxb is the total number of LMXBs in 
the galaxy, our model would predict that Lx, total /L op t 
varies linearl y with Sn- While not the same relation- 
ship found bv I White et al.1 l|2002|) . we find that the data 
are equally well fitted by our model. It has also been 
shown that A^lmxb and Lx, total in early- type galaxies 
are proportional to the stellar mass of the gala xies with 
a wea k dependence on the morphological type J^Jilfanov 
12004^1 . Our model would also predict that -/Vlmxb and 
Lx, total are proportional to the stellar mass of the galax- 
ies with some dependence on the GC-specific frequency 
that could yield the mor phological type d ifference seen. 
It should be noted that iGilfanovl l)2004f) suggest that 
statistics and inaccurate calibration of the near-infrared 
mass-to-light ratio may cause some of the morphological 
dependence seen in their data. Our model would predict 
that the dependence on Sn would be of the order of a 
factor of 1.5 given the galaxies in the l^ilfanov, ( 2004) 
sample, similar to what was measured. However, the in- 
clusion of late-type galaxies may actually decrease the 
significance of this effect given the expected difference in 
the field LMXB population formation. A further study 
encompassing early-type galaxies with a larger range of 
Sn values is warranted. 

The high (Lx > 5 x 10 37 erg s _1 ) observed luminosi- 
ties of the LMXBs in early-type galaxies has been used 
to suggest that such systems could not have formed from 
primordial binaries since the X-ray lifetimes of such sys- 
tems would be much shorter than the stellar age of the 
galaxy. The general argument is that at 5 x 10 37 erg s _1 
and with a 1 Mq companion, the system would have an 
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X-ray active phase of less than m 100 Myr. There are 
many difficulties with this conclusion. First, the lag time 
between binary formation and LMXB activity can be 
quite long. If the companion is a 1 M© donor that starts 
mass transfer after the donor has turned off the main 
sequence, then the binary was formed more than 10 Gyr 
ago. For systems that start mass transfer on the main se- 
quence, the time between the formation of the binary and 
the start of mass transfer is dictated by either gravita- 
tional radiation or magnetic braking, both of which have 
chara cteristic timescale s of the order of a few gigayears 
( fi .g.. lPfaM et^l20f)l . 

It is also important to recognize that the mass-transfer 
rates observationally inferred may not accurately reflect 
the long-term mass-transfer rates due to short-timescale 
fluctuations, such as the disk ionization instability or 
X-ray irradiation of the co mpanion (e.g., Ivan Paradiisl 
1199a iBuning fe Ritterl l2004) . As a result of such effects, 
intrinsically low (~ 10~ 10 M Q yr _1 ) mass-transfer rate 
systems can masquerade as high mass-transfer rate sys- 
tems. All of the BH LMXBs in the Milky Way, as well 
as a comparable number of NS systems, exhibit X-ray 
outburst phases in which the maximum luminosity is a 
significant fraction of the Eddington limit for the com- 
pact object. These systems, commonly referred to as soft 
X-ray transients (SXTs), are thought to have low mass- 
transfer rates from the donor such that the accretion disk 
is uns table to the disk ionization instability ijvan Paradiisl 
1996). In SXTs, the systems undergo outbursts lasting 
days to months (and even years in some cases) and then 
return to a quiescent state. The recurrence time can be 
as short as a year or longer than the history of X-ray as- 
tronomy. While the disk instability model works well in 
predicting which systems are likely to undergo this phe- 
nomena, it provides no clues as to the length of the X-ray 
outbursts or recurrence times. Other transient phenom- 
ena, such the effect of X-r ay irradiation on the donor 
(e.g. JBuning fe Ritterll200^1 ." can also cause intrinsically 
low mass-transfer systems to have observed luminosities 
much higher than expected. 

3. GLOBULAR CLUSTER CONTRIBUTION TO FIELD 
LMXBS 

It has been suggested that GCs may contribute to the 
field LMXB population throug h ejection of LM XBs or 
disruption of the clusters fe.g-. lWhite et alJ2002f) . These 
alternatives have been invoked for a number of reasons 
including the expected short X-ray active lifetimes of 
LMXBs formed from primordial binaries and the simi- 
larities in the X-ray properties of the field and GC popu- 
lations of LMXBs. But it is important to recognize that 
once LMXBs leave the GC environments, through either 
ejection or disruption, they will evolve like primordial 
binaries and would have similar X-ray active lifetimes. 
Therefore, GC formation of field LMXB scenarios require 
a constant input of LMXBs into the field to account for 
the observed population. 

The same dynamical processes that form LMXBs in 
GCs would be responsible for the ejection of these sys- 
tems. Therefore, the number of LMXBs ejected from 
GCs in any galaxy would be proportional to the number 
of GCs or S N M*. If all of the field LMXBs were due to 
ejections, we would expect no dependence of /lmxb,gc 
on Sn- Our data rule out this scenario at greater than 



the 99.99% confidence level. While not required, we can 
not rule out that some small fraction of the field LMXB 
population was ejected from GCs with the current data. 
The upper limit on the fraction of the field population 
due to ejections is dependent on Sn- For Sn — 2.0, we 
find that < 20% of the field population of LMXBs could 
have been ejected from GCs. 

Theoretical studies of GC disruption find that a signif- 
icant fraction of the GCs in a galaxy will be disrupted 
with the exact fraction s howing some dependence on 
the g alaxy properties fe.g.. lVesperinill2000llFall fc Zhand 
2001). If we assume that disruption of a GC does not 
change the probability for LMXB formation in the clus- 
ter, then to fully account for the number of LMXBs seen 
in the field of early-type galaxies, the number of dis- 
rupted clusters must be between 0.5 and 4 times the num- 
ber of GCs currently observed. Such numbers push the 
limits of w hat is expected from current models for GC 
disruption ijVesperinil 12000) . but LMXBs from disrup- 
tions cannot be ruled out. One difficulty however, is that 
a significant fraction of the disruptions may occur early 
in the history of the galaxy, in the first < 3 Gyr. At early 
times, the dynamical interaction s in GCs may be less effi - 
cient at producing LMXBs fe.g.. lDavies fc Hansenll998|l . 

4. DISCUSSION 

We have shown that the observed relationship between 
the fraction of LMXBs found in GCs and the GC-specific 
frequency in early-type galaxies is consistent with the 
field population of LMXBs forming in situ via primor- 
dial binary formation. To test our model, better data for 
both /lmxb.gc and Sn are required. In current studies, 
only a fraction of the LMXB population in early-type 
galaxies can be compared to the GC population because 
of differences in the field of views of the instruments. 
More thorough coverage of the Chandra field of view by 
GC studies will improve the statistics for /lmxb.gc- The 
values for Sn are also subject to observational uncertain- 
ties because of uncertainties in the estimates of both the 
number of GCs and the optical magnitude of the galax- 
ies. Of particu lar note is the rec ent revision in Sn for 
NGC 1399 fsee lDirsch et^l200l . 

A measurement of a time-dependence in the number of 
field LMXBs per unit galaxy mass would also verify our 
model. This measurement requires good estimates of the 
stellar age of galaxies. In addition, understanding how 
both stellar age measurements and LMXB populations 
are affected by complex star-f ormation hist o ries is neces- 
sary to interpret the results. iTrager et al.l l)2000|) found 
that galactic age estimates from single stellar population 
models are heavily weighted by young stars. No com- 
parable study of the theoretical prediction for LMXB 
populations has be en performed. As noted previously, 
IWhite et all l)2002j) found that the total X -ray luminos- 
ity attributed to LMXBs in early-type galaxies was not 
correlated with stellar age, but their analysis included the 
contributions from both field and GC LMXBs. The GC 
LMXBs are not expected to be dependent on the stellar 
age of the galaxy, and therefore the inability to separate 
GC and field LMXBs, due to the low spatial resolution 
of AS CA, may explain wh y no age dependence was found 
in the IWhite et all lj2002|) study. 

The spatial distribution of GC s is generally more ex - 
tended than the galaxy light (e.g,. IAshman fc Zepi l998). 
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TABLE 1 

Globular Cluster and LMXB Connection Data 



Galaxy 


Sn 


Nlmxb.gc 


J^LMXB a /LMXB.GC (%) 


Refs 


NGC 1332 


2.2±0.7 


9 


30 


30t^ 


1.2 


NGC 1399 


5.1±1.2 


26 


38 


68±« 


3,4 


NGC 1553 


1.5±0.5 


2 


11 




5.6 


NGC 3115 


1.6±0.4 


9 


36 


25+1 


7,8 


NGC 4365 


4.3±0.6 


18 


37 


«3 


5,9 


NGC 4472 


3.6±0.6 


30 


72 


421^ 


10,11 


NGC 4486 (M87) 


14±4 


60 


96 




12,13 


NGC 4649 


6.7±1.4 


22 


47 




5,14 


NGC 4697 


2.1±0.3 


32 


89 


36±5 


15,16 



a The total number of X-ray point sources covered by the field of view of the 
optical globular cluster catalogs. 

References: (1 ) iKundu fc Whitmord <200lT), (2) iHumphrev Buotd <200l) . 
(3) IDirsch et alj 120031) (4) lAngelini etaT] 120011) . (5) IKissler-Patisl 119971) . 
(61 ISarazin et alj 120031) . (71 IPuzia et alj I2004ir ~(8l IKundu et al.l 120031). "(9 ) 
ISivakoff et alj 120031) (101 IRhode & Zepfl 120041). f l 11 IMaccarone et alj 120031. 
(12) IMcLaughlin et at] 119941) . (13) IJordan etaU 120041) . (14) IRandall et alj 
120041) . (151 Jordan et al. 2005, in prep., (16) Sivakoff et al. 2005, in prep. 



In addition, GC spatial profiles flatten at small radii. 
Our model predicts that given the different formation 
sites for GC and field LMXBs, there would be similar 
radial distribution differenc es seen in the LMXB popu- 
lation. IKundu et al.l ((2002) found that the radial pro- 
files of field and GC LMXBs were roughly similar, al- 
though their study was limited in the number of sys- 
tems (30 and 42 for GC and field LMXBs, respec- 
tively) and the radial extent sampled. Better cover- 
age of the LMXB and GC populations is necessary, es- 
pecially at the inner and outer regions of the galaxies 
where the GC and optical light profiles are the most 
different. Two additional effects must also be consid- 
ered: i) how supernova kicks will affect the distribution 
of primordial binaries in the field, and ii) how the pref- 
erence of GC LMXBs to be located in red, rather than 
blue. GCs llAngelini et all200lHKundu et all200l 12001 
iSarazin et all 120031: IJordan et alJ 120041) affects the GC 
LMXB distribution. 

Given the different formation mechanisms for field and 
GC LMXBs, it is reasonable to assume that the distri- 
butions of the physical parameters (i.e., orbital period, 
donor mass, and mass transfer rates) will be different. 
This, in turn, could yield differences in the observed lu- 
minosity function. Yet how different the luminosity func- 
tion will be is unknown. At present, no theoretical study 
of the predicted physical parameters of GC LMXBs ex- 
ists. Such a study is essential in determining if a lumi- 
nosity function difference is expected and/or measurable. 



Also, a better understanding of how transient behavior 
affects the observed luminosity function is also necessary. 
Multi-epoch observations of LMXBs in early-type galax- 
ies will provide more information on the transient behav- 
ior of LMXBs, although longer observations are required 
to distinguish between real transient behavior (luminos- 
ity changes of more than 3 orders of magnitude) and the 
factor of ~ 10 variabilit y commonly seen in persistent 
LMXBs (|Lin et alJl2001l and references therein). 

Finally, it has been suggested that the lack of spec- 
tral differences between field and GC LMXBs supports 
a com mon origin for these systems l|Maccarone et all 
20031). There is no observational evidence of differences 
in the spectral pa rameters of field and GC L MXBs in the 
Milky Way (e.g.. iChristian fc Swanklll997|) . Therefore, 
it is unlikely that any spectral differences would be seen, 
particularly given the low spectral sensitivity of studies 
of LMXBs in early- type galaxies. 
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